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Abstract. Let M be a smooth manifold and S a spray defined on the convex 
cone C of the tangent bundle TM. It is proved that the only non-trivial fc-jet 
approximation of the exact geodesic deviation equation of <S, linear on the 
deviation functions and invariant under arbitrary local coordinate transforma- 
tions corresponds to the Jacobi equation. However, if linearity in the deviation 
functions is not required, there are differential equations whose solutions admit 
fc-jet approximations and are invariant under arbitrary coordinate transforma- 
tions. As an example of higher order geodesic deviation equations we study 
the first and second order jet geodesic deviation equations for a Finsler spray. 



1. Introduction 

Given a connection on a manifold M and two geodesies of the connection, it is 
a fundamental problem in mathematical physics to determine the relative motion 
between the pair of geodesies. This problem is of great interest in theories of 
gravity and in astrophysics. The solution to this problem for Lorentzian manifolds is 
provided by the solutions of the Jacobi equation, the Jacobi fields along Lorentzian 
geodesies pi [24]. 



Although the geodesic equation is not linear, the Jacobi equation is a linear differ- 
ential equation. Therefore, the solutions of the Jacobi equation can only describe 
the deviation between the geodesies under some assumptions, that sometimes are 
not applicable to interesting physical situations. This fact motivates the search for 
consistent generalizations of the Jacobi equation. Currently, there are two frame- 
works for generalizations of the geodesic equation extensively used in astrophysics 
and general relativity: 

• The theory developed by Hodgkinson and Ciufolini and based on 
the linear rapid deviation hypothesis and whose fundamental equation we 
denote as the generalized Jacobi equation (see equation (jlip below). 

• The formalism developed by Bazahski and others, based on equations for 
higher order jet fields (see for instance P fT^)- 

In this work we consider the problem of the general covariance of the first approach 
in the general setting of geodesies of arbitrary sprays. The methodology used 
involves a generalization of the second approach to arbitrary sprays. 

One difficulty of the first approach is the identification of the geometric character 
of the solutions of the generalized Jacobi equation. It was proven that the solutions 
cannot be tensorial [9]. Our suggestion in this paper is to consider /c-jet fields for 
the solutions of the generalized Jacobi equation. This is a natural choice, because 
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any local, natural, smooth differential operator can be represented as a /c-jet field 
(see Peetre's theorem, see [R,, P- 176). As a result of our analysis, we conclude 
that for an arbitrary spray, if the solutions of the generalized geodesic deviation 
equation are approximated by /c-jet fields, the generalized Jacobi equation from 
Hodgkinson 12 and [7 is not general covariant. The proof of this fact is based on 
the analysis of the hypothesis under which the equation is obtained, showing that 
under some specific coordinate transformations, such hypothesis breaks down. 

The structure of this work is the following. In section 2, we introduce some technical 
definitions that will be used through the paper. In particular, we define the notions 
of general covariance, approximation scheme, exact geodesic deviation equation, 
generalized geodesic deviation equation and jet approximations. We state the main 
problem considered in this paper. Our results from the analysis of the problem 
is presented in the form of Theorem A and Corollary B. In section 3, we develop 
further the notion of jet approximation to deviation equations. Section 4 contains 
the proofs of the main results and also Lemma H.W In section 5, we discuss higher 
order geodesic deviation equations. The idea is applied to connections associated 
with a Finsler spray. We make some digression on the consequences that the higher 
order deviation equations have for Finslerian cosmology and Finslerian relativity. 
Finally, in section 6, we discuss the results and their relation with previous work. 



2. Preliminaries and results 

Let M be an n-dimensional smooth manifold and C C TM a sub-bundle of codi- 
mension zero such that each fiber is convex cone of T^M. Given a coordinate 
chart {U,{z'^ : U M", n = l,...,n}), coordinate indices are indicated by Greek 
characters, and run from 1 to n = dim{M). Repeated up and down Greek indices 
indicates the sum over 1 to n if nothing else is stated. We will assume that all the 
geometric structures required in our considerations are smooth. 

Geodesic equation of a spray. A second order differential equation (or semi- 
spray) is a smooth vector field S S TTC such that tt^S{u) = u, \/u G C, with 
IT : TC — ;> C being the canonical projection. By using a local frame, a semi-spray 
can be expressed in the form 

(1) '^(^'2^) = ^'a|r-'^'(^'^)a^- 

If the semi-spray is homogeneous, its components are homogeneous functions of 
degree 1 on the y coordinates and therefore S'^{x,y) — y'^^^-g^^ ='■ u'^ ^ai^^v)- 
Then iS is a spray. In this case, the integral curves of the spray S are invariant 
under affine reparameterizations. 

The geodesies of S are the projection to M of the integral curves of S. Also, the 
homogeneity condition for the spray S implies 

n{x,y)^ yP^^^^^, [x,y)€T^M, f,,,, ^ 1, ...,n. 

Using this relation, the geodesic equation of S can be written as 

(2) z'^ + r^(z,i)i'^ = 0, fi=l,...,n. 



The functions r(^(z,i:) are the non- linear connection coefficients (see for instance 

m)- 
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There is associated to 5 a canonical, linear, torsion-free connection V of the tangent 
bundle ttc : TC ~> C (see for example |20j . Chapter 1). The non-zero formal 
Christojfel symbols are defined by the expressior0 

Then the geodesic equation ([2]) is written as 

(4) + r(;,(z,i)i''i^-o, ^l = l,...,n. 

They define an afhne connection if they live directly on M, that is, if each coefhcient 
Ti^pix) e -F°°(M,R), for each fi,iy,p^ l,...,n. 

We consider deviation equations associated with geodesic equations only, and we 
will assume that the connection is symmetric, that is, the relation 

(5) ^^.a=^L, P,'^,<J^l,...,n 

holds in any coordinate system. However, the method of analysis can be extended 
to more general differential equations. 



Geodesic deviation equations. Let us consider two nearby geodesies x: I ^ M 
and X : I ^ M, / C M and assume, in order to simplify the treatment, that the 
images x{I) C M and X{I) C M of both geodesies are on the chart domain U C M . 
Let : / — > K, ^ = 1, n} be the coordinate displacement between the geodesies 
defined as 

(6) ^'^(s) := a;^(s) -X^(s), sel. 

Since x : I ^ M and X : I ^ M are solutions to the geodesic equation, we have 
the relation 

(7) 't + K,{x + ^ + [x'' + r) + r ) - KAX) X- X- = 0. 

This relation is referred to as the exact geodesic deviation equation. One can think 
the relation ([T]) as a constraint. Let us introduce the function F : M'^" — > M", 

(8) 

(^M^^M^^M) + r(;^(x + e,x + c) (x'' + c") + r ) - T^,AX) x" x"" . 

Then the exact geodesic deviation equation ([7]) reads as F(^(s), ^(s), ^(s)) = 0, 
where F is interpreted as a function along the geodesic X : I ^ M . 



Notion of general covariance. A pseudogroup of transformations of a topologi- 
cal space T is a collection of local homeomorphisms of T such that the composition 
of transformations is compatible with the collection (for a formal definition see for 
instance [TS], p. 1 or [53] )■ A pseudo sub-group F' of a pseudogroup F is a subset 
of F which is a pseudogroup itself. Pseudogroups are of relevance to the problem 
that we are considering because the notion of general covariance of a differential 
equation can be formulated as the invariance of the differential equation under the 
action of certain pseudogroup of transformations. 



Note that some of the Christoffel symbols are chosen to be zero for any natural frame. In par- 
ticular the associated with covariant derivatives along vertical directions, for instance V a Z = 

for any vector field Z G TTC. 
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Definition 2.1. Given a smooth function G : M"i R"" with m e N, G is 
compatible with the pseudogroup r°°(M"i) of local diffeomorphisms $ : — 
of class C°° if it is invariant by the action of the pseudogroup, i.e., 

Go$=G, V$e r°°(M"i). 



It is assumed that the variables on which the function G depends hold a representa- 
tion of r°°(IR"i). For the case that we are investigating, the relevant pseudogroup 
is r°°(R"). We consider the class of C°° coordinate transformations since we will 
require Taylor's expansions of functions on AI of arbitrary order. 

Approximations of the exact deviation equation. The relation (O is com- 
patible with the action of each element of the pseudogroup r°°(M"). However, the 
equation ([7]) is non-local in M, since the definition of {^'^ : / M, fx = 1, ...,n} 
involves for each s two points on the manifold M . This implies a geometric charac- 
ter which is difficult to understand. One way to avoid this difficulty is by allowing 
the following approximations, 

• The Christoffel symbols TP,^{X + ^, X -f ^) are approximated by Taylor's 
series in terms of the functions and 

• Assuming that are infinitesimal, i.e, the monomials {^'^^'', 

etc...} are dropped out in a given equation (in the case of geodesic 
deviations, the equation is = 0, with F being the function ([S])). 

The relation ([7]) under the above approximations yields to the following expression, 
that can be thought as an ordinary differential equation 

OTP- ■ ■ 

(9) Cr + -Q^{X) e.X'^ X" + 2 VP,{X) i\X^ = 0. 

This is a local equation in the sense that only data defined along the central geodesic 
A : / — > is necessary in its formulation. Also, equation ^ is linear and the set 
of solutions defines a finite rank vector bundle of Jacobi fields along A : J — >■ M. 
Equation ([9]) is a non-explicit covariant way of writing the Jacobi equation jl8] . It 
is compatible with the pseudogroup r°°(R") and one can write it in a explicitly 
covariant way as 

(10) V^V^,J + i?(J,A)A = 0, 

where J(s) = ^ii^)l^\x{s) ^ '^x{s)M and R{J, A) is the Riemann type curvature 
endomorphism of the connection V determined by J, A S Tx(s)^^- Then under 
the above hypothesis, the solutions of the exact deviation equation ([7]) can be 
approximated by the solutions of the covariant differential equation (fTU)) . that are 
sections of the jet bundle Jg(M, A/) along the central geodesic X : I M. 

If one makes the assumption that only the deviation functions {S,^ , /z = 1, n} are 
infinitesimal (that is, that only the monomials j^''^'^, /x, = 1, n} are negligible), 
the approximation of the relation ([7]) for any V connection yields the qeneralized 
Jacobi equation [12J, 

(11) kP + r^, {ikpx" + kpk"'^ + ^^k" (a" + /i") [x" + /i") = 0. 

Being a non-linear equation, the unknowns k^" have an obscure geometric inter- 
pretation. This paper is motivated by the problem of understanding the geometric 
nature of pTjl and in particular the issue of the general covariance for geodesic equa- 
tions associated with arbitrary smooth sprays for general sprays. This includes the 
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problem for affine and Finsler sprays. The covariance problem is equivalent to the 
compatibility of the equation pT|) with the action of the pseudogroup r°°(]R"'). In 
practice we will investigate a weaker notion of general covariance, based on the 
invariance under the action of certain pseudo sub-groups of r°°(]R"). 

Jet approximations. Motivated by the difficulties in determining the geometric 
character of (fTTj) . it is reasonable to assume a specific geometric character for the 
solutions of the differential equation (|lip . The fields {k^} cannot be tensorial, since 
for instance the space of solutions of equation dTTI) cannot be a vector space (for 
a detailed discussion of the non-tensoriality of the fields {k^} see [Sj)- Also, the 
solutions are not represented by irreducible representations of the Poincare group as 
the fields described in [25] . Other natural candidates for the solutions of equation 
(fTTj) are sections along the curve X : / — )• M of a fc-jet bundle tt^ : Jg (/, M) — > M. 
These are smooth maps \1/ : / -> Jq (/, M) such that the diagram 

(12) JS{I,M) 




commutes. Some relevant fields defined along X : I ^ AI are of this type. The most 
relevant example is the case when k — 1 and Jq{I, M) corresponds to the tangent 
bundle tti : TM M and the sections {S,*^ : I ^ J^{I, M)} along X : / -)> M are 
the Jacobi vector fields by J(s) = ^i{s)i^\x{s), solutions of the Jacobi equation 
®. 

In order to introduce fc-jet fields along the geodesic AT : / — s- A/ it is convenient 
to work with 1-parameter geodesic variations associated with the pair of geodesies 
X,x : I ^ M. This is a map 

A : (-60, eo) x I ^ M, 

such that for each e the function A(e, •):/—> Af is a geodesic and such that 
A(0, s) — X{s) and there is an e G (— eo,eo) such that A(e, s) = x{s). Then the 
local coordinate functions A''(e, s) are expanded in the first variable e G (— eojCo) 
by Taylor's theorem. As a result, we find two applications for our problem. First, 
it allows us to approximate the fields {k^ : / ^ M, ji = 1, ...,n} in terms of k- 
jet sections and perform an analysis order by order in e of the covariance of the 
equation ((TT|) . Second, it allows us to define fc-jet sections along X : I ^ M 
and a perturbative approach to the problem of finding deviation equations. This 
corresponds to a generalization of Bazahski's theory (see for instance [5] or |14)). 

By embedding the geodesies X : I ^ M and a; : / Af in the ribbon A((— eg, eo) x 
/) C M we can approximate the functions , fi — 1, ...,n}, solutions of the exact 
deviation condition ([7]), by using Taylor's expansion 

i=0 

for a fixed e G (— eoi ^o)- Indeed, it is useful for a variational interpretation of the 
deviation functions ^'^ to consider the following ribbon coordinate Junctions, 

(13) : (-eo,eo) x/^K, (e, s) A^(e, s) - A''(0, s). 



The error in the approximation 

S^(e,s)^l]^(e,s) 
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is given by the remainder term of the approximation, which is of order 6*^+^ . There- 
fore, the error in the approximation £,^{s) by ^^^(e, s) is of order 6*^+^. 

The approximation — > Q'^ is an application of Taylor's theorem. One problem 
considered in this work is if the further approximation 

(14) K(s),^-l,...,n}^{r!^(e,s)}, 

where {k^(s), /j, = is a solution of the equation ([TT|) and e S (— eojEo) is 

fixed as before, is consistent with the pseudogroup r°°(M"). For arbitrary sprays, 
the answer to the problem is negative. This is a question more related with the 
approximations involved in obtaining the differential equation (jlip than with the 
structure of the equation itself. Indeed, the non-covariance of the approximation is 
generic for deviation equations (note that the equations are not necessarily geodesic 
equations) obtained under the same hypothesis than the generalized Jacobi equation 

m- 

Let us assume that in the 1-parameter family of geodesies A : (— eo,eo) x / — > M, 
the parameter e € (— eo, ^o) is invariant by the action of each element of r°°(K"). 
Although sometimes such parameter can be constructed explicitly (for instance in 
the case of Riemannian manifolds e can be the distance between the geodesies), 
in general we need only to assume that it exists for the geodesic families that we 
consider. 

Approximation schemes. Let us consider a geodesic variation A : (— eo, eo) x / — >■ 

M associated to the original pair of geodesies x : I ^ M and AT : / — > M by an 
embedding (see Proposition^^. 

Definition 2.2. The monomial ^^{s)^'^{s) is negligible at order k if for the Taylor 
expansions up to order k the monomial ri{^(e, s)Sl^(e, s) is of order 0{e^^^). 

Similarly, the definition applies to monomials of type C'^C, £,^5,'^ : to and 
other monomials. The property of being negligible is covariant, since we have 
required that e is an invariant parameter. 

The different approximations found in the literature to the exact geodesic deviation 
equation was the motivation for the following 

Definition 2.3. An approximation scheme is a set of negligible monomials at order 
k of the free algebra generated by the monomials {^^, i^", iJi,v, p = 1, ri}. 

Example 2.4. The following three examples are considered in this work: 

• Trivial approximation scheme, where any of the monomials 

is not negligible. 

• Linear approximation scheme, where all the monomials 

are negligible. 

• Linear rapid deviation scheme, where from {£^5,^ , C^^'^, C^'^, /i, = 1, n} 
the only negligible monomials are , p,^ — 1, ...,n}. 

We will consider the approximations of the type (jl4p and then apply a particular 
approximation scheme. Note that for fixed e, the field r2^(e, •) is determined by 
a A;-jet along X : I ^ M. The full scheme of approximations that one makes is 
compiled as follows: 
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(i) The deviation functions {^^ : / — > M, = 1, n} are substituted by the 
ribbon coordinate functions {S'^ : (— eo,eo) x / — > K, = l,...,n}. The 
justification of this fact is based on Provosition 13 . 1 1 of section 3. 
(ii) The ribbon coordinate functions {S^ : (— eo,eo) x / — >■ M} are approx- 
imated by the Taylor expansions {il^ : (— eo,eo) x / — > M}. This ap- 
proximation can be done with arbitrary accuracy for an appropriate fc, by 
Taylor's theorem. 

(Hi) One assumes a particular approximation scheme. In particular the linear 
rapid deviation approximation scheme is equivalent to the approximation 
-j> fi^, = 1, n in the relation ([7]). 

These three assumptions together with the embedding proved in Proposition 13.11 
imply that the fc-jet fields {£7^ : (— eo, eo) x / M} are also an approximation up 
to order e^'^^ for the solutions {k^ : / M, /i = 1, of the equation (ITT]) . 

Method and results. The consistency criterion to be checked against general 
covariance is that the error in each of the approximations introduced must be 
bounded by or be of the same order than the error in the approximation {^^ 
0^(e, s), /i — l,...,n}. Note that the value e has been fixed. However, in order 
to have a variational interpretation, we consider the family defined by the variable 
parameter e G (— eo,eo)- If e < e, the error will be also bounded. Therefore, we 
consider errors bounded or of the same order than powers of e. This reformulation 
allows to speak of order / in e and make our criteria equivalent to the following: 
for fc-jet approximations, all the error must be of order fc -|- 1 in e, which is the 
error in the approximation S'^(e, s) by ri^(e, s). The error in considering 
negligible must be comparable to the error of considering fi^ri^ negligible. To 
be compatible with the Taylor's approximation by $1^, this error must be of 
order e*'"''^, i^'^^k — C'(e*'+^). Similarly, the following conditions must hold for the 
monomials: 

(15) n^,n'^ c,0{e'^+'), n^.n-^, 0{e''+'), n^^n-^ 0{e'^+'), etc... 

The conditions can be checked order by order in A: G N. As a consequence 
of such analysis, we show that using fc-jet field approximations to the solutions of 
equation (fTT|) one can test the compatibility with r°°(]R") of the approximation 
(fT5)) of the exact Jacobi equation. In particular, we have the following 

Theorem A. Let 5 G F TC be a spray. If the solutions of the associated gener- 
alized deviation equation ([TT|) are fc-jet fields along the geodesic X : I ^ M, then 
either: 

• Equation pT|) is not compatible with the pseudogroup r°°(IR"), or 

• The separation between the geodesies x,X : I ^ M is not small. 

From the proof of Theorem A in section 3, it follows that 

Corollary B. Let S G FTC be a spray. Given two geodesies X : I ^ M and 
X : / — > i\f , the only fc-jet approximation to the exact deviation equation ([7]) linear 
on K and compatible with F°°(M") is for fc = 1 and corresponds to the standard 
Jacobi equation. 

The strategy to prove these results is to consider some specific non-empty pseudo 
sub-group r°°(A, fc) C F°°(IR") and check which approximation schemes are com- 
patible with F°°(A, fc) for each fc > 0. The pseudogroup r°°(A, fc) is more manage- 
able than the full pseudogroup F°°(K"). The negative outcome of our check implies 
the results. 
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3. Jet approximations of geodesic deviations 



Let us consider the spray S and the associated connection V. The following result 
shows that locally, two nearby enough geodesies x : I ^ M and X : / — ?> M can 
be described as a members of a 1-parameter geodesic variation for a short time 
interval J C /, with / = [0, a]. 

Proposition 3.1. Let x : I ^ M and X : I ^ M be two geodesies such that 
x{0) and X(0) are connected by an extensible transverse curve 7. Then there is a 
1-parameter geodesic variation A : (— eq, eq) x / — > Af such that the central geodesic 
is A(0, s) — X{s) and A(e, s) — x{s) for some e £ (— eo, eo). 

Proof. The two initial points x(0) and X{0) can be joined by the connecting curve 
7 : [0, e] — > M with 7(0) = X{0) and 7(e) = a;(0). Indeed, one needs to extend a bit 
the curve 7 to include [0, e] in an open interval (0 — u, eo), with e < cq. To obtain 
the desired geodesic variation we construct appropriate initial conditions along 7. 
First, the tangent vector X{0) is parallel transported along 7, defining a vector 
field along 7 denoted by ^'(e). A similar operation can be done for i(0) but along 
the inverted curve —7 from x{0) to X(0), generating a vector field x'{e). Let us 
consider the linear combination of vector fields along 7, 

Z(e,0) = i (ex'(e) + {e~e)X'{e)) € T^^,)M, 1^0. 

By Picard-Lindelof 's theorem, the set of initial values (7(e), Z{e, 0)) determines an 
unique geodesic A(e, ■) : [0, Smax{<^)] C I ^ M for some maximal / 3 Smaxi^-) > 0. 
It is clear that 

I 3 Smax mzn{sniax(e),e e [0, e]} > 0, 

by compactness. By continuity, the same is true for a small enough and Smax 
defined as 

Smax := r7im{s,„ax(e), e e (0 - ct, e - a)}. 

Therefore, we have constructed a geodesic variation A : (0 — ct, e+a) x [0, Smax] — > M 
with A(0, s) = X{s), A(e, s) = x{s) for s G [0, Smax]- By a convenient reparameter- 
ization of 7, the parameter in the variation can be redefined in the interval (— eo, eo) 
and still keep A(0, s) — X{s) and A(e, s) — x{s) in the new parameterization of A 
with the required properties. □ 

It is not essential which parallel transport we use in the construction of the vector 
field Z. One can use for instance the parallel transport of the connection V on 

the total lift of X{e) to TC, and then push-forward the corresponding vector field 
along 7. Another option is to use the parallel transport of any Riemannian metric 
defined on M. 

Corollary 3.2. Let x,X : I ^ M and A : (— eo, eo) x I ^ M as in Provosition \S.l\ 
Then there exists a e G (— eo,eo) such that the relation £,^{s) = A''(e, s) — A^(0,s) 
holds. 



We can use Taylor's expansions on e of the smooth maps A^ : (— eo, eo) x / ^> M to 
approximate C'^(s). In particular, by application of Taylor's theorem up to order k 
to the function A^ it follows that 

l_ 9J A^(e,g) 
fc! del e=a 



(16) A^{e,s)^J2 



3=0 



a'=+iA^(e,; 



(fc + 1)! de^+^ 
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with e(s) £ (0, |e|). Taking derivatives respect to the parameter s in (|16p one obtains 

^ '^ ^ fc! del e=o (fc + l)!V ae'^+i e=.-(.)A 

We define the foUowing fields r2^(e, s) and £'^'^^(s) by the relations 
(17) 

^ ^' i—l 

For the fixed value of e such that x{s) — A(e, s). The corresponding time derivative 
are 

J7^(e,s):= A'^(6,s)-A''(0,.)- -^^^^+i(s)= ^ |er(.), /i=l,...,n. 



1=1 



The fields {Cl'(s) := — a j''^ , fJ- = j = l,...,fc} live along the central 

" e=0 

geodesic A(0, s) = X{s). £7^ are fc-jet bundles along X : I M. On the other 



hand, the remainder fields 

(18) e^+i(s) ^^Ffi 



e=£(s) 



do not live along the geodesic X : I M. Motivated by this fact, one can define 
the following set of transformations. 

Definition 3.3. Fixed a 1-parameter geodesic variation A : (— eo,eo) M, r°°(A, fc) 
is the collection of diffeomorphisms (p : M" — )■ M" such that under the action of 
r°°(A, fc) the functions (S'^, il^) transforms to (S^,f2^) with the constraint 

(19) S^(e, s) - n^,{e, s) = S^(e, s) - n^,{e, s),Vfi = 1, n. 

Proposition 3.4. r°°(A, fc) is a non-empty pseudo sub-group of r°°(R") with the 
product composition of mappings. 

Proof. To prove that it is non-empty, let us a consider transformation $ : M" — s> R" 
in the same fc-jet than the identity transformation Id : R" ^> M", but with different 
^-jet for / > fc + 1. Then the restriction of $ to A((— eq, sq) x I) c U C M provides 
an example </? G r°°(A, fc). The fact that it is a pseudo sub-group it is proved 
by checking that the axioms of pseudogroup hold [IS]. That r°°(A, fc) is a pseudo 
sub-group of r°°(R") is direct. □ 

Given the above Taylor's expansions of order fc and a particular approximation 
scheme, a second order differential equation described by an algebraic expression 
G(^,^,^) — can he approximated by another algebraic expression -D(^,C,0 = 
by equating to zero in G(^, ^, ^) — the monomials that are negligible. Such 
approximation is consistent with fc-jet expansions if the error G — >■ 13 is bounded or 
of the same order than the approximation {^'^j {Vl'^{e, s)}: if the error G — ?• I? is 
not bounded by the error {^^} — ?> {i^^}, then one cannot say that is a solution 
to G(f2^, f2^, f^^) =0, which is the fc-jet approximation to deviation equations that 
we are interested.. This imposes a constraint on the schemes of approximations 
that are compatible with the pseudogroup r°°(K"). Therefore, the problem of the 
compatibility of a approximation scheme with r°°(R") is translated to check if 
the error in the approximation of the monomials ^^^"^ by r2^(e, s)n^(e, s) is 
compatible with r°°(R"). In practical terms, it is enough to check the compatibility 
with {r°°(A, fc),Vfc G N}, for all the monomials that are negligible. 
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4. Proofs 

In this section we prove Theorem A and Corollary B. We assume that ah the 
functions and expressions are smooth. We start proving the following Lemma, 

Lemma 4.1. For any general covariant functional expression G{^,^,^) = 0, the 
only approximation schemes compatible with r°°(R") are: 

(i) The trivial approximation scheme, where none of the monomials 

{^'^r, eC", , er^", •••} are negligible, 
(ii) For Taylor expansions with k — I and in the linear approximation scheme, 
where all the monomials {S,'^^", C^C"^, S.^^", ■} are negligible. 

Proof. The trivial approximation scheme is clearly compatible with r°°(R"'). Let us 
consider the linear approximation scheme. For k=l, the monomials {^^S,'^ , ^.^S," , 
^M^'^^^M^'^^P) can be approximated by monomials n'^n'(,n';;n'(,h'^n\nP) 
with an error of order without imposing any restriction on the character of 
the functions il^(e, s). Therefore, such approximation scheme is compatible with 
r°°(]R"). Indeed, for fc = 1 the monomial ^^^^ is negligible as the following short 
calculation shows, 

n^^ie, sMie, s) = e2er(s)er (s), Ve £ (-eo, eo). 

This implies that the term ^i{e, s)ili{e, s) is of the same order in e than the error 
in the approximation S''(e, s) — > ni{e,s) in any coordinate system, for each e G 
(— eo, eo)- Since we assume that the remainder is negligible, one can neglect the term 
V.i{e, s)ili{e, s). Similarly, for the monomial il^l{e, s)Cl'[{e, s) one has the relation 

n^,ie,sme,s)^ e(l^{s)eCA^)=^'^n^)iU^) = 0{^). 

If we require that for k=l the monomial (^^C" is not negligible, then from the above 
relation it follows that ^^{e, s)tli{e, s) must be of order i. Thus, it must be a 
smooth function C^'^ : / — > R such that 

(20) ^'i{e,s)tll(e,s) = C^''{s)-^ 

holds. However, this requirement is not compatible with r°°(A, 1). A simple con- 
stant scaling in local coordinates — Az^ (that induces an element of r°°(A, 1) 
by smooth restriction to the band A((— eo, e) x /)) implies a change = A^^ (the 
fields {^f, ^ = l,...,n} are tensorial, since they are the components of a Jacobi 

vector field). If we take the constant value A = e, it follows that also = A^f and 
therefore 

fi5'(e, s)h'i{e, s) ^ €^n''^{e, s)tl'i{e, s) = a"'{s)e^ ^ = C^"'{s)e. 

This is a contradiction with the relation (P(7| , except 

(i) If C^''(s) — 0, implying that the Jacobi field gfjr is zero. In this case, the 
condition (j20l) is covariant, but the associated curvature endomorphism TZ 
must be zero, or 

(ii) The parameter e = 1, which is a contradiction with the general require- 
ment that fJ.,!^ = l,...,n} are negligible. To see this, let us consider 
e to be the distance function between the curves x(s) and X{x), for a 
given Riemannian metric on M . Then it is clear that if e = 1 there are 
geometric configurations where ^^^^ (for instance) can be of order 1, and 
therefore negligible iff is negligible. 
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A similar argument follows for other monomials. 

For fc > 2, the argument is analogous. Let us consider the Taylor approximations 
of order k for and the corresponding expansion of the monomial 



1 



A; 



\_ , 9JA^(e,s) 



^ IT 



9^A''(e,s) 



kl 



5e' 



Each of the monomials must be negligible, which implies that they must be of the 
same order e^^+i than the remainder g'^_^i{e,s), 

^U^f. c)'A^(e,s) 

c=0 



dei e=o 9e' 



9eJ 



for some positive p. For j = I = \, this statement is not covariant, since implies 
the condition 



^ > V de e=o de 



aA''(e,s) aA''(e,s) 



de 



= 0(ef) 



for a positive integer p. In a similar way as for the case k — 1, one can prove that the 
condition (PT|) is not covariant under the action of the pseudo sub-group r°°(A, k). 
Therefore, non-trivial approximations schemes such that ^^^"^ are negligible only 
could work for fc = 1 or one of the following possibilities hold: 

(i) If e = 1, which is a contradiction with the fact e is small, 
(ii) If all the fields {C^} are null. This possibility is impossible for arbitrary 
initial conditions on the geodesies. 



□ 



Proof of Theorem A. For fc = 1 it follows from Lemma 14.11 that the approxi- 
mation scheme where all the monomials {^''^''} are negligible is compatible with 
r°°(A,fc), but then the rest of the monomials S^^S,", etc... must also be negligi- 
ble. Thus, if there are no restrictions on the curvature endomorphisms, for fc = 1 
the assumptions under which equation ([TT]) is approximated from the relation ([7]) 
does not hold in arbitrary coordinate systems and therefore, the approximation 
scheme leading to the generalized Jacobi equation is not invariant under the action 
of r°°(A, fc). Therefore, for fc = 1, the only approximation scheme compatible with 
r°°(A, fc) is the linear approximation scheme, leading to the Jacobi equation ©, 
that we know is compatible with the pseudo sub-group r°°(A, fc). For fc > 2, it 
follows from Lemma HTT] that the approximation scheme where ^'^^'^ is negligible is 
not compatible with r°°(A, fc) except in the situations when e is not small (which 
is a contradiction with the requirement that {^^S,"^ , = 1, n} are negligible. □ 

Corollary B follows directly from Theorem A. 

5. Non-linear approximation schemes and Finsler space-time higher 
order geodesic deviation equations 

We have analyzed the possible approximation schemes where are infinitesimal. 
However, Theorem A leaves open the possibility for the existence of alternative ap- 
proximation schemes where the monomials {'^^^'', /i, = 1, n} are not negligible. 
In this case, one can still use Taylor's expansion on the parameter e and fc-jet field 
approximations. The method is a direct generalization of Bazahski's theory from 
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Riemannian setting to general sprays. If 5 is a spray defined in a cone C ^ TAf , 
let us consider the expansions (|17p and its derivatives and insert them in the exact 
deviation equation ([7]) . The connection coefficients are also develop in the variable 
e, obtaining a finite polynomial in e. Then the relation ([7]) can be written as a 
formal series, 

oo 

fe=i 

Equating to zero each term, one obtains a hierarchy of ordinary differential equa- 
tions, 

(22) Gfc(S^S^S'^) = 0, k = 1,2,3,... 

The equation obtained from the first order Gi(S^, S'^, S^) = is the Jacobi equa- 
tion of the connection associated with .S. Higher order deviation equations are 
obtained by equating to zero the expressions G/c(S'', S'', S'^) = for k — 2,3,.... 
These higher order geodesic deviation equations are still valid even if S^^S^" are not 
negligible. 



Higher order geodesic deviations in Finsler geometry. The spray S must 
not necessarily be affine. In particular, one can consider Finsler sprays living on the 
slit tangent bundle N — TM \ {0}. Then we can adopt the following generalization 
of the notion from J. Beem of Finsler space-time to arbitrary signature metrics [3]: 

Definition 5.1. A pseudo-Finsler structure on the slit tangent bundle N is an 
homogeneous of degree 2 smooth map L : iV — > M such that the vertical hessian is 
not degenerate on N. 

Definition 5.2. A spray S E TTN is a Finsler spray if there is a Finsler function 
L : N ^ M. such that the integral curves of S are the total lift of the geodesic curves 
of L. 



Given a Finsler space-time, there is an Ehresmann connection (see for instance 
determined by L, defining a decomposition 



(23) 



TN = 'H(BV, 



where V = ker^di:). An adapted frame to the horizontal- vertical decomposition is 
determined by the smooth tangent basis for T„iV for each u G N: 

(24) 

where iV* j are called the non-linear connection coefficients. Given a tangent vector 
X S TxM and u G Tr^^{x), there is an unique horizontal tangent vector h{X) € 
TuN with dTT{h{X)) = X (horizontal lift of AT). 

The pull-back bundle Tr*TM is the maximal subset of the cartesian product N x TM 
such that the diagram 



TT*TM ■ 



■TM 



N ■ 



M 
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commutes. There is a Chern's type connection i^llj defined on ■k*TM and it can 
be shown that there are three types of curvatures associated with the spray, the 
Riemannian type or /i/i-curvature, the /iw- vertical curvature and the vertical or vv- 
curvature. The Riemann type curvature is defined to be the tensor field R with 
components 

(25) 

Kp^x, y) = {^Kp - -^Kp - K>^L - Kx^'-pu) y). ^.P.^ = 1, n. 
The /iw-curvature of the spray S is the tensor whose components are 

(26) Pl^p^ix,y) -.^ -^{x,y), ^i,v, p,a = 1, ...,n. 

By the condition of homogeneity of the function L, one can see that the /iv-curvature 
vanishes along a geodesic, — 0. In the case that the associated covariant 

derivative is zero along vertical directions, the vertical curvature is automatically 
zero. This is the case for the Chern connection, whose only non-trivial curvatures 
are the Riemann type curvatures and /i?;-curvature. Note that although the tensors 
R and P are linked by Bianchi identities [T], in general R does not determine the 
/iw-curvature tensor P. 

In the following, we present the first and second order geodesic deviation equations 
for Finsler space-times. First, let us introduce some notation. The covariant vector 
J2 is defined by the expression in components J2 = Ji'^^ where the components 
are [H] 

(27) J!^ = e2 + Kp^iil- 

The vertical lift of a tangent vector Z e T^M to T(^^_y-jN is denoted by = Z>^^. 
We denote by Cz the Lie derivative along Z acting on tensor fields. 

Proposition 5.3. Let F : N ^ ^ he a Finsler function and consider the expansion 

Then 

• The first order geodesic deviation equation of a Finsler spray is the Jacobi 
equation, 

(28) V^V^ei+ 0. 

• The second order geodesic deviation equation is the non-linear differential 
equation 

v^y^j2 + R^{J2,x)x = Ve,i?(x,6)^-V;^i?(x,6)Ci + 4i?(x,a)(v^a) 

(29) -I- C^^^^,yP{{\/^^^r,X)X + C^^.^,y.R{Ci,X)X. 



Proof. Both equations can be obtained simultaneously by developing equation ([7]) 
in powers of e. The first equation is a generalization of the Jacobi equation for 
Finsler space-times [llj . obtained by grouping together all the terms proportional 
to e in the expansion of the exact geodesic deviation equation ([7]). When this is 
done, one obtains the condition 

dV^ ... ... 

+ ^(^' X) cf^" X" + 2 nAx, X) x-^er = o. 

A re-arrangement of this expression one obtain equation ([25)1 . Equation follows 
from the equality G2 = in front of the term in the exact deviation equation. 
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Regrouping terms that arc proportional to e , one obtains the expression 

OTP ... ... . . OTP ■ ■ ■ 

+ -Q^ix, X) c^x" X- + 2 T^.ix, X) x^ii + 2 vp^i\i{ + 4^(x, X) cfx'' r 

9^r^ . . . . 9^r^ . . . 9^r'^ . . 

Re-arranging the terms, the above expression is equivalent to 
(30) 

V;,V;^ J2 + i?^(j2,x)x= V5,i?(x,a)^-V;^i?(x,ei)ei + 4i?(x,a)(Vj^6) 

. . a^r*' . . . a^r^ . . 

The first line corresponds to the second order deviation equation in the case of an 
affine connection and is covariant (see [M]). The second line is related with the 
/iti-curvature and is intrinsically a non-affine contribution. The functions {^^ , /i = 
1, n} do not define the components of a vector field along X : I M. In order 
to define an associated vector field, one can consider the covariant derivative Vj^^i 
and the functions 

(32) i^x^iT - er + rp{x, x),pCiXp. 

In order to re- write the second line in expression pip in a covariant way, one can use 
normal coordinate systems, where the r(^^(a;, y) = at a fixed poinilj. In normal 
coordinates and using the definition of the /iw-curvature, the second line can be 
written tensorially as 

-f£(^,).i?(ei,i:)i:. 

This is a tensorial relation and equation ([25]) holds in any coordinate system. □ 

Riemann-flat Finsler space-times. In Finsler geometry with positive definite 
Finsler metrics, it is a difficult problem to find metrics of non-Riemannian or 
Minkowski metrics with Riemannian curvature tensor zero R — 0, but the /in- 
curvature is non-trivial (see pQ; pg: 328. The existence of singular examples of 
such spaces for the Chern connection follows from a result from Zhou [21]). In 
Lorentzian signature, one way to relax the problem is to assume that the domain 
where the metric L is smooth is not the slit tangent bundle N and allow singular 
regions where L is not smooth or is degenerate. A Finsler space-time with R = 
and P ^ will be called Riemann-flat Finsler space-time. 

In Riemann-flat space-times the first deviation equation (|28p reduces to 

(33) V^V^Ci=0, 
and the second deviation equation ([29]) reduces to 

(34) V^V^J2= q^,yPiiS/^^,r,X)X. 

Even if i? = 0, the equation Vj^Vjj^^i = is not equivalent to = 0, since in 
principle there is not a coordinate system where in a whole open neighborhood 
TP,p{X,X) = holds. Equations (|33p and ([M]) have an interesting interpretation 
in Finslerian cosmology (see for instance |17)'). since induces an additional bending 



•^For the existence of such normal coordinate systems, see for instance Chapter 5. Note 
that at the origin, the exponential map of the Chern's connection is only . However, this is 
enough for our argument. 
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for geodesic curves. This is caused by the local anisotropy of the Finsler space-time 
and is encoded in both fields and J2. Also, it is interesting to remark that if 
ordinary matter is determined by the curvature R, the curvature P must be linked 
with the vacuum structure of the space-time. 

The local time inversion operation exchanges the sign of the derivatives respect to 
the inversion of the parameter t: for a function / G J- {I, M), the local time inversion 
operation maps ^/ ~^/- Note that in order to define this operation it is not 
necessary to speak of a global time orientation. The effect of the action of the 
time inversion operator on the geodesic X : {a, b) — > M is to reverse the geodesic 
TX : {a,b) M,X{b~t + a) — X{s). The induced action on the tangent vector of 
X is ^ Thus, if is a solution of the equation (l33l). ^^(b-t + a) is 

also a solution of the reversed differential equation ([55]) along the reversed geodesic. 
It follows easily that for Riemann-flat space-times, the first order geodesic deviation 
equation is invariant, 

Proposition 5.4. For Riemann-flat spaces, the Jacohi equation p3p is invariant 
under the action of the local time inversion operation. 

The second order geodesic deviation equation of a generic Finsler space-time is 
not invariant under 7". For instance, this is the case if the space-time metric is 
non-reversible. Therefore, the tensor field 

(35) V -.TTM xTTMTTM ^TTM, {X,Y, Z) ^ C^y^yy Pii"^ xY)'" , X)Z, 

when evaluated on the geodesic {X, X) is sensitive to the time inversion. This time 
asymmetry in the geodesic deviation can be interpreted as a geometric source of 
irreversible evolution. 

Another striking consequence of equations and ([M]) is related with the valid- 
ity of the Einstein equivalence principle in Finsler geometry. To show this, let us 
consider the Lorentzian case. In Lorentzian geometry, there are smooth normal 
coordinates such that the Levi-Civita coefficients r(^p(a;) are zero at a given point 
p G M and the coordinate functions are smooth. By a convenient choice of the co- 
ordinate systems, one can disregard the effects of the curvature term i? in a small 
enough space-time region, and the geodesic equation describes the dynamics of a 
free of force, free-falling point particle. This construction is justified by the fact 
that normal coordinates are smooth and by the fact that in Lorentzian geometry, at 
the selected point p G M, the connection coefficients r(^p(x) can be approximate as 
r^p(x) « in a small neighborhood, by continuity. This follows by the possibility 
of approximating the curvature i? w in a small enough region, again by conti- 
nuity. Note that this possibility is taken by guaranteed in the formulation of the 
week equivalence principle. However, the equations p3p and p4p forbid in general 
the existence of smooth normal neighborhoods. In particular, normal coordinates 
are only of class C^. Thus, the above argument by continuity of the approximation 
ri^p{x) w in a small open neighborhood is not valid for Finsler geometry. The 
physical consequence of this fact is that Einstein's equivalence principle does not 
hold for generic Finsler space-times and in particular for non-trivial Riemann-flat 
space-times. The obstruction that avoids the existence of smooth normal neigh- 
borhoods is that the tensor P is non-zero. Indeed, the existence of normal smooth 
neighborhoods implies that P — 0. In the positive case, those spaces are called 
Berwald spaces. The analogous concept of Berwald space for Finsler space-times is 
straightforward. 

Definition 5.5. A Finsler space-time (M, L) is of Berwald type iff for each p £ M 
there is a coordinate system such that the connection coefficients y^pip) vanish. 



16 ON K-JET FIELD APPROXIMATIONS OF GEODESIC DEVIATION EQUATIONS 

We have indeed the foUowing 

Proposition 5.6. Let (M,L) be a Finsler space-time. Then if Einstein's equiva- 
lence principle holds, then the Finsler space-time (M, L) is of Berwald type. 

Proof. If Einstein's equivalence principle holds, it is possible to define a coordi- 
nate system where the effects of gravitation are locally eliminated (a freely falling 
laboratory coordinate system). This is equivalent to the fact that the connection 
coefficients of the Chern type connection are smooth and zero at one given point. 
Therefore, the space must be of Berwald type. Conversely, if the Finsler space-time 
is of Berwald type, Einstein's equivalence principle holds automatically. □ 



6. Discussion 

In the Lorentzian case, equation has been applied extensively in astrophysics 
(see for example [THl IH O H] and references there) and in astrometry (see for ex- 
ample [S]). The affine case was investigated in [21 . Although (TTlT) looks as being 
non-covariant, it was argued that the generalized Jacobi equation (jlip is equivalent 
to an explicitly covariant ordinary differential equation. Indeed, it was written in 
terms of the covariant derivative of the Levi-Civita connection V and its curvature 
tensor R (sec for instance, [HI [311]) in a manifestly covariant way. However, the 
known arguments for such equivalence between the covariant and the non-covariant 
form of the generalized Jacobi equation (|lip contain hypotheses on the curvature 
R or other hypotheses whose general covariance character is quite doubtable to 
be valid. For instance, in |T2] there is a relevant passage concerning the order of 
magnitude of the approximations which specially significant, in particular when 
some second order terms in are disregarded, in the partial linear approximation 
scheme ( see [T^], pg 365-366). This is in contradiction with Lemma (j4.ip 'l. since 
the hypotheses that second order terms are negligible is consistent with general 
covariance only for k — 1 and for the linear approximation scheme. In particular, 
for rapid deviation schemes it is not longer true that the square of the deviation 
functions are negligible. In [i] the authors have made the assumption that the 
components of the Riemannian curvature tensor are small. Without more specifi- 
cation, this is a non-covariant statement. In [7], the author has used the assumption 
that certain parallel transport between the respective points x{s) and X{s) of the 
geodesies does not depend on the curve connecting them. However, this is only 
true under restricted curvature conditions (in particular, for trivial holonomy con- 
nections). In view of the above criticism, we should conclude that the equivalence 
between the equation (|lip and the covariant versions found in the literature are not 
valid for general Lorentzian space-times and for general connection sprays. Thus, 
the applicability of equation ([TT|) as approximation to ^ is restricted to special 
coordinate systems or for geometries with special constraints in the curvature such 
that the holonomy of V restricted to the chart U is trivial. 

It was argued in 22J that there is not a consistent generalization of the geodesic 
deviation in the rapidly deviation regime. Schutz's analysis relies on the prescrip- 
tion that the geodesic curves are joined by geodesies and is restricted to Riemann 
normal coordinates. In contrast, we have only required an initial curve connecting 
the initial points a;(0) and ^(0). This method extends the conclusions of [22] to 
more general pairs of geodesies, not necessarily with image on normal coordinate 
domains. Moreover, we have seen that there are general covariant generalizations 
of the Jacobi equation (see for instance [21 [H] and the discussion the deviation 
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equations for Finsler sprays). However, one needs to abandon the requirement that 
the monomials {£,^£,'^, ^ly — 1, ...,n} are neghgible. 

The apparent contradiction between the conclusions of this paper and [22 with the 
use of the generalized Jacobi equation in astrophysical applications can be dis- 
entangled if the generalized Jacobi equation (1111) is not general covariant. Indeed, it 
is enough to realize that it is not compatible with the pseudo sub-group r°°(A, k). 
If equation (jlip is not general covariant, the arguments that have been given in [2 2) 
in normal coordinates are not necessarily applicable in Fermi coordinates, that are 
the coordinates where the applications of equation ([TT|) have been developed (see 
for instance [3J [5J [5J [21] ) • Note that the geometry of Riemann normal coordinate 
neighborhoods is very different than the geometry of Fermi coordinate neighbor- 
hoods, and that the second type of geometry (roughly speaking, a cigar open set) 
is not as sensitive to detect rapid deviated orbits as the spherical normal coordi- 
nate geometry (roughly speaking, an isotropic ball type open set). Therefore, by 
the restriction that the geodesies stay in a tubular neighborhood around the cen- 
tral geodesic, the linear rapidly deviation scheme leading to (jlip is consistent with 
the chosen geometry. Then it is applicable only if the geodesies do not deviate 
too much, which is in contradiction with the motivation of the equation. On the 
contrary, the rapid deviation approximation scheme is not consistent in a generic 
normal coordinate system, since the geodesies can deviate in an appreciable and 
fast way. This argument explains why the calculation in [22 is performed in normal 
coordinates and also the limitation of the use of the generalized Jacobi equation to 
Fermi coordinate calculations. 

Equation (jlip is covariant under affine coordinate transformations. Since the trans- 
formation between two Fermi coordinate systems are affine transformations, the 
generalized Jacobi equation (jll[) is covariant under local coordinate transforma- 
tions from Fermi to Fermi coordinates [9 . However, as an approximation to the 
exact deviation equation ((T]) , the generalized Jacobi equation (ITT|) fails to be general 
covariant, since the hypothesis of the approximation scheme break down in arbi- 
trary coordinates. This is consistent with our analysis, since the pseudo sub-group 
of local afHne coordinate transformations is in general different than r°°(A, k). 

There are alternative frameworks to the generalized geodesic deviation for going 
beyond linearization in the analysis of geodesic deviations. Bazahski's theory [2], 
is a convenient framework to investigate geodesic deviations beyond the Jacobi 
equation. In such formalism, one can formulate an hierarchy of general covariant 
differential equations for the fields {^j , fi — 1, n, j — 1, fc}. Bazahski's theory 
was used extensively in the investigation of geodesic motion in general relativistic 
space-times (see for instance in [T3] and in subsequent works of those authors). 
We have found that such methods can be extended to connections determined by 
arbitrary sprays S G FTC 

We have applied the generalized Bazahski's theory in the case of a Finsler spray, 
obtaining the classical Jacobi equation in Finsler geometry and the second 
order deviation equation (1291) . Such higher order geodesic deviation equations can 
be used if the deviation between geodesies are not negligible. In the Finslerian case 
and for the Chern type connection, not only the Riemannian curvature is involved, 
but also the /iw-curvature tensor P is significant. Therefore, one expects that the 
/i?;-curvature could play a relevant role on determining the topology of the pseudo- 
Finsler spaces. For positive Finsler spaces, this is of relevance for the Berwald- 
Landsberg problem |10| and have implications for Finsler space-time cosmologies. 
We also have seen that generically in Finsler space-times the Einstein equivalence 
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principle and time inversion symmetry does not necessarily hold. These effects have 
been demonstrated for Riemann-flat Finsler space-times, but one expects similar 
effects for general Finsler space-times with R ^ 0. 
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